EXPLICIT DIFFERENTIAL CHARACTERIZATION OF 
PDE SYSTEMS POINTWISE EQUIVALENT TO Y xiixi2 = 

1< h,h < n>2. 

JOEL MERKER 

ABSTRACT. In Lie 1 883 , as an early result, Sophus Lie established that a second order 
ordinary differential equation y xx = F(x,y,y x ) is equivalent, through an invertible 
point transformation (x, y) i— » (X(x, y), Y(x, y)), to the free particle equation Yxx — 
if and only if the right member F is a degree three polynomial in y x , namely there exist 
four functions G, H, L and M of (x, y) such that F can be written as 

F(x, y,y x ) = G(x, y)+y x - H(x, y) + {y x f ■ L(x, y) + {y x f ■ M (x, y), 

and furthermore, the four functions G, H, L and M satisfy two second order partial 
differential equations: 

= -2G yy + -H xy - ^L xx + 2(G L)y - 2G X M - AGM X + ^HL X - ^HH y , 

= ~\ H vv + \ L *v ~ 2M *x + 2GM V + 4G y M - 2{HM) X - ^H y L + ^LL X 

In IM2004a) , this theorem was generalized to systems of Newtonian particles with m > 2 
degree of freedom, i.e. with one independent variable x and m > 2 dependent variables 
(y 1 1 V 2 1 ■ ■ ■ i y" 1 )- ln m i s P a P er . we generalize S. Lie's theorem to the case of several 
independent variables (x 1 , x 2 . . . , x n ), n > 2, and one dependent variable y. Strikingly, 
as in IM2004al . the (complicated) differential system which corresponds to the above 
two second order partial differential equations is of first order . By means of computer 
programming, this phenomenon was discovered in IBN2002I . IN2003I in the case n = 1. 
In IBi2003l . IM2003I . the general case n > 2 was handled. 
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§1. Introduction 

This paper, a direct continuation of |M2004a|, provides a summarized proof of the 
following statement, labelled as Theorem 1.23 in the Introduction of |M2004a|. All our 
functions are assumed to be analytic. 

Theorem 1.1. (n = 2: IBN2002I . IN2003I : n > 2: IBi2003l . IM2003n Let K = E 

or C, let n G N, suppose n > 2 and consider a system of completely integrable partial 
differential equations in n independent variables x = (x , . . . ,x n ) S K" and in one 
dependent variable y G K of the form: 

(1-2) y x i,. x to{x) = F n ' n (x,y(x),y x x(x), . . . ,y x n(x)) , j x ,h = 1,.. .n, 
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where F J1,J2 = F J2 ' J1 . Under a local change of coordinates (x,y) t— > (X, Y) = 
(X(x, y), Y(x, y)), this system (1.2) is equivalent to the simplest system Y X j 1X j 2 = 0, 
j\,j2 = 1, ■ ■ • ,n, if and only if there exist arbitrary functions Gj 1 ,j 2 , i, L k * and 
M kl of the variables (x, y), for 1 < Ji, J*2, ki < n, satisfying the two symmetry condi- 
tions Gj lt j 2 = Gj 2 ^ 1 and - 2 = Hj*h> suc ^ l ^ at eauat ' on (1-2) is of the specific 
cubic polynomial form: 
(1.3) 

n ( 1 1 \ 

y x31x32 = G juh + Vx«i ( + 2 Vxn L jl + T, Vx*» L jl + Vx .11, 2 ) » 

fci=i 



forji,j 2 = l,...,n. 

We refer the reader to |M2004a| for an extensive introduction and for a more sub- 
stantial bibliography. Applying E. Cartan's equivalence algorithm (see IG 19891 and 
IOL 19951 for modern expositions), the general case n > 2 of the above theorem has also 
been established in IHal937l . where the construction of a projective connection associ- 
ated to a second order ordinary differential equation achieved in ICa 19241 was extended to 
several variables. Theorem 1.1 was re-discovered in |BN2002|, in |N2003|, in |Bi2003| 
and in |M2003|, thanks to partial parametric computations of the Hachtroudi-Chern ten- 
sors in n > 2 variables, which were described in a non-parametric way in ICh 19751 
(see §1.8 below). 

It may seem quite paradoxical and counter-intuitive (or even false?) that every sys- 
tem (1.3), for arbitrary choices of functions Gj lt j 2 , if? 1 - 2 , L k ^ and M kl , is automati- 
cally equivalent to Y X j 1X j 2 = 0. However, a strong hidden assumption holds: that of 
complete integrability. Shortly, this crucial condition amounts to say that 

(1.4) D xi3 (F^)=D xj2 (FM>), 

for all ji,j2,j3 = 1, • • • j n, where, for j = 1, . . . , n, the D x , are the total differentiation 
operators defined by 



W /; . : 7777 ' 7T7, • X 1 



9x3 d y i^i dVx 



These conditions are non-void precisely when n > 2. More concretely, writing out (1.4) 
when the F J1,J2 are of the specific cubic polynomial form (1.3), after some nontrivial 
manual computation, we obtain the complicated cubic differential polynomial in the vari- 
ables y x k labelled as equation (1.25) in lM2004a| . Equating to zero all the coefficients 
of this cubic polynomial, we obtain four familes (F), (IF), (IIF) and (IV) of first order 
partial differential equations satisfied by Gj 1J2 , H^^, and M kl : 



(F) < - G juj2iX j 3 - G jl!hiX j 2 + Y H ji,h G ki,h - X! H h,h ( '' 
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(IF) 



r o = 57 g 



(HI') 4 



<5f G 



if 



fei 



2 ^31 ,33 ^ 2 ~ 2 33 + 



1 1 

2 X! G k>,33 L j2 - 2 S h Zl G k2M L ): 
k 2 = l k 2 = l 

1 n n 



fc 2 =i 



fe 2 =i 



Errfcl irt 2 _ rrfci rrfc 2 

^2,33 31,32 n k 2 ,j 2 n jl,h- 



k 2 = l 



k 2 =l 



= 



o-ee 2 



„fe 2) Tjk<,(ri _ r fe <r(2) TT k a{l) 

33 h,32,y h 31,33,2/ 



, 1 r fc <r(a) r fe <j(i) ^ r fc ^(2) r^Ci) , 

' 2 / • ii _ 2 ? ' 31 + 

,1 rfc^(2) r*-V(l) 1 c fc ^(2) r fc <T(l) , 

+ 2 Jl 32 _ g 1 



■n. 



k 3 = l 



fcs=l 



1 1 ' 

2 d 3i 2w W fe3,33 - 2 « ^ fc ^' 2 33 

fcs=l fc 3 = l 



1 k c 

- 2 5 32 



ka=l 



2 

1 A 



ir k °{2) jk 3 _ ^ c^(i) 
Z-^ k 3 ,j 3 3i 2 ' 



Err fe <r(2) ,-fcg 
W fc 3 ,32 ^3i 



fc 3 =l 



y H k 3 A«_lA(D y £ 

2 33 / / 31,32 k 3 2 32 / i 31,33 

fc3=l fe3=l 



fe "(2) 

fc 3 



(IV) < 



CT6S3 



1 rkcr( 3 ) ,fccr(2) j^k^(l) ^^^(3) ,&cr(2) j^kcr(l) 



2 33, 3i 32, a 2 J2 ' -? 1 



33, V 



_|_^fco-(3)>fc<r(2) J^ k a{l ) _ ^(3)^(2) J^ k "W 



32, 31 



'32 , 31 



2:^3 



33, 31 



X-J2 



/ I K4,33 33, 31 / < K4,32 

fe 4 =l 

r 

fco-(l) ,feo-(3) 



fe 4 = l 



1 cfc aC i),k„ C s) 
4 31 , 33 



E r fe^(2> r fe 4 1 s-k CT (i) ,fc CT (3) r^* 2 ' r fc i 
^fc4 ^32 4 °3"l, 32 ^4 

fc 4 =l 



fc 4 =l 



These systems (F), (IF), (IIF) and (IV) should be distinguished from the systems (I), 
(II), (III) and (IV) of Theorem 1.7 in IM2004a| , although they are quite similar. Here, 
the indices 31,32,3%, k\, fca, £3 vary in {1, 2, . . . , n}. By 62 and by S3, we denote the 
permutation group of {1, 2} and of {1, 2, 3}. To facilitate hand- and Latex-writing, partial 
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derivatives are denoted as indices after a comma; for instance, Gj 1 j 2iX j 3 is an abreviation 
for dG jul2 /dx h . To deduce (F), (II'), (HI') and (IV) from equation (1.25) of |M2004al . 
it suffices to observe that every cubic polynomial equation of the form 

n n n 

o = a+ y B ki- y x *i + Y Y c **te • y* ki y* k2+ 

ki=l ki=l k 2 = l 

(1.0) 

v 7 n n n 

+ Y D kiM-k 3 ■ Vx k i Vx k 2 Vx k z 

k\ — l k 2 — l k-2 — 1 

(as for instance (1.25) in (M2004a| ) is equivalent to the annihilation of the following 
symmetric sums of its coefficients: 

= A, 
= B kll 

= C*i,k 2 + Cfc 2 ,fei I 

= Dk 1 ,k 2 ,k 3 + Dk 3 ,k!,ka + B>k 2 ,k 3 ,k 1 + Dk 2 ,k 1: k 3 + Dk 3 MM + Dk 1 ,k 3 ,k 2 - 

for all k±, k 2 , k 3 = 1, . . . ,n. 

In conclusion, the functions Gj lt j 3 , H^* , and M fel in the statement of Theo- 
rem 1.1 are far from being arbitrary: they satisfy the complicated system of first order 
partial differential equations (I'), (IF), (HI') and (IV) above. 

Our proof of Theorem 1.1 is similar to the one provided in IM2004al . in the case of 
systems of Newtonian particles, so that most steps of the proof will be summarized. In the 
end of this paper, we will delineate a complicated system of second order partial differ- 
ential equations satisfied by Gj lt j 2 , H£ L k ^ and M kl which is the exact analog of the 
system described in the abstract. The main technical part of the proof of Theorem 1 . 1 will 
be to establish that this second order system is a consequence, by linear combinations and 
by differentiations, of the first order system (T), (IF), (IIF) and (IV). Before proceeding 
further, let us describe rapidly a second proof of Theorem 1.1, confirming its validity. 

1.8. Confirmation of Theorem 1.1 by the method of equivalence. Let us summarize 
the str ategy of IBN2002I . IN2003I . IBi2003l . In ICh 19751 . inspired by M. Hachtroudi's 
thesis |Hal937 ], S.-S. Chern conducted the equivalence algorithm associated to the sys- 
tem (1.2). His computations are essentially equivalent to M. Hachtroudi's (though with 
much less information), but they are done in a non-parametric way. It was deeply known 
to Elie Cartan and it is now well known by modern followers (cf. for instance IG 19891 . 
IGTW 19891 . IHK1989I . IFe 19951 . IOL 19951 ) that achieving parametric computations in 
an equivalence problem is an extremely hard task, often devoted to computers nowadays. 
In fact, the so-called Cartan Lemma (|OL1995 1, p. 26) was essentially devised by E. Car- 
tan as a shortcut, alongside his groundbreaking progresses, after some experiences of 
explicit hand computations. The trick in this lemma is to bypass the (often too long) para- 
metric computations, nevertheless keeping track of some important informations. Nowa- 
days, to study an equivalence problem, one often achieves the non-parametric computa- 
tions by hand, leaving the explorations of parametric computations to a computer. In this 
respect, the thesis IN2003I of Sylvain Neut is extremely interesting, since it implements 
the general equivalence algorithm as a Maple package. 

However, computer programs achieving formal algebrico-differential computations 
with a general number n of variables are still undevelopped. Consequently, it is in- 
teresting to achieve the computations of ICh 19751 in a parametric way, because at the 
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= 0, 1 < Ji, 32 < n > 2 



end, the vanishing of all the invariant tensors that one obtains in the final {e}-structure 
would show explicitely under which condition the system (1.2) is equivalent to the sys- 
tem Y X j lX j 2 = 0. To our knowledge, this problem is considered as open in the field 
of symmetries of Cauchy-Riemann submanifolds of C'\ a subfield of the mathematics 
area called Severa] Complex Variables. However, no specialist of the subfield seems to 
be aware of the old reference !Hal937l . a text only read, apparently, by S.-S. Chern, who 
was a student of E. Cartan contemporary to M. Hachtroudi. 

Let us expose briefly how one obtains the final {e}-structure attached to the equiv- 
alence problem associated to the system (1.2). Since we have not been able to follow 
everything in the ancient text IHal937l . we follow ICh 19751 with the only change that 
we do not introduce the imaginary number i = ^/— T in our structure equations and also, 
we make some minor changes of sign; in I Chi 9751 . the computations are conducted over 
K = C with the idea that they apply to Complex Analysis, but they do hold as well over 
K = R, or even over an arbitary commutative field of characteristic zero equipped with a 
valuation (in order to provide K-analytic functions). 

To begin with, consider the following family of 2n + 1 initial differential forms: 



lo := dy - y x n dx 13 , 



(1.9) 



LO 
LOn 



dx c 



:=dy x .-J2 Fa,P d^- 







Here, the Greek indices a, (3, 7, S, e, £, etc., run from 1 to n. Sums J2a=i are abbreviated 
as J2 a - I n me paragraph preceding the statement of Theorem 1.7 in IM2004al . we explain 
why we cannot use coherently the Einstein summation convention throughout the paper. 
For the same reason, we shall abandon this convention in the present paper. However, 
we would like to mention that in (1.9) above and in (1.10), (1.11), (1.12), (1.13), (1.14) 
and (1.15) below, the summation convention applies coherently, so that the reader may 
drop the sums if (s)he is used to. 

Following ICh 19751 . define the local Lie group consisting of (2n + 1) x (2n + 1) 
matrices of the form 



(1.10) 







.10 



where u a and u a are n x 1 vectors close to the zero vector, where u % is a n x n invertible 
matrix close to the identity matrix, where u is a scalar close to 1 and where u'^ denotes 
the inverse matrix of Up. With this group, define the initial G -structure, consisting of the 
following collection of (2n + 1) differential forms, called the lifted coframe: 



(1.11) 



( lo := u- LO, 



to a := u a ■ lo ■ 



6 



JOEL MERKER 



These forms depend both on the "horizontal" variables (x l , y) and on the "vertical" 
("group", "fiber") variables u, u a , u a and Ug. The introduction of this lifted coframe 
may be justified as follows. 

It is the very first step of the method of equivalence to check that there exists a trans- 
formation [x 1 , y) i — ► (x l ,y) of the completely integrable system y x a x a = F a 'P in the 

coordinates (x 1 , y) to another completely integrable system y x a x n = F ' in other coor- 
dinates (x 1 , y) if and only if there exist functions u, u a , u a , up* and vp* of (x l , y), which 
depend on the functions (x l ,y), on their partial derivatives with respect to (x l ,y), on the 
y x i and on the F 1,s , such that the following matrix identity holds: 




where the (2n + 1) differential forms (uj, cu , uj a J are defined similarly as in (1.9), in the 

barred coordinates. Of course, it is understood that in the left hand side of this matrix iden- 
tity, the variables (x 1 : y) are replaced by their values with respect to the variables (x l , y) 
(in the language of modern differential geometry, one usually speaks of "pull-back"). By 
a more careful examination of the explicit expressions of the functions u, u a ,u a , up* and 



Va in terms of the F 1 ' 6 , one observes that v@ = uu' a ( IM2003| ). Thus, the collection of 
differential forms (u>, u> a , U) a ) is defined modulo multiplication by a matrix of functions 
of(x t ,y) which is of the specific form (1.10). Based on this preliminary observation, 
the general procedure of the equivalence method ( IG1989I . IOL 19951 ) associates to the 
system (1.2) the lifted coframe (1.11). 

Applying the exterior differential operator d to u>, to u> a and to u> a and absorbing the 
torsion, it is shown in ICh 19751 that the initial structure equations may be written under 
the form 



(1.13) 



du = oj a A uj a + co A if, 

a 


du a = if^Aujp+uJa Atf + u A<p a , 





where if, ifp*, f a and f a are modified Maurer-Cartan forms. We notice that there are no 
torsion coefficient to normalize. In fact, the choice of v% := u u'^ made in advance above 
corresponds to having achieved a first normalization implicitely. 

Since the dimension of the Lie symmetry group of the system (1.2) is always finite and 
in fact bounded by n 2 + An + 3 ( IHal937l . ICM 19741 . ICh 19751 . ISu2001l , IGM2003I ), 
according to the general procedure of the method of equivalence (|G1989|, |OL1995|), 
one simply has to prolong the initial G-structure. One could also apply the Elie Cartan 
involutivity test to deduce that it is necessary to prolong. 

Now, we summarize the remainder of ICh 19751 very rapidly. After one prolongation, 
two normalizations and one supplementary prolongation, the final {e}-structure is of the 
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following form: 

doj — ^2 uj°~ A u a + w A (f, 

ot 

duj a = uj 13 A ipp + uj A ip a , 



duj a — fa A OJp + W a A ip + uj A ip a , 



dip = \^ W " A If a — ^ WaA^"+wA!/l, 

d ^ = J2 J2 s pp • A u " + J2 R p-> • u A ^ + £ T r ■ w A 



+ w Q A + Sp ^2 ^ A + ui/3 A ip a + ^2 f'p A V™ + \ $g ' u A ip, 



(1.14) { 



■ W A uj. 



+ i^Q^^A^+^L^ 



LJ A UJ/3 + 



(3 7 



+ y> A ip a - <P0 A <p + i <^ Q A ip, 



d<p a = Rl ■ UJ P A uj 1 + - Qi ■ uj A top + P c 



p ■ u) A w + 
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dtp = ^ • uj a A ujp + ^2 H <x ' w A u a + ^ A"" • uj A uj a + 

a j3 a a. 

+ ip Aip + 2^2 ip a A <p a . 



These structure equations incorporate 8 families uj, oj a , uj a , ip, ipg, p a , tp a , ip of differ- 
ential forms, of total cardinality n 2 + An + 3, together with 8 invariant tensors S^p, Rp~, 
Tp 1 , Q<£, L a P, P a0 , H a and K a having some specific index symmetries that we shall 
not use. 

Applying the exterior differential operator d to these 8 families of structure equa- 
tions (1.14), one verifies that the seven invariant tensors Rp 7 , T^ 1 , Q^, L a/3 , P a p, H a 
and K a are in fact functionally dependent on the fundamental tensors Sp°, namely they 
are certain coframe derivatives of the Sg". In the (very similar) context of the equivalence 
problem associated with a Levi non-degenerate hypersurface of C™ +1 , this computation 
was achieved by S.M. Webster in the Appendix of ICM 19741 : in the precise context of 
the equivalence problem associated to the system (1.2), this computation is not achieved 
in »Ch 19751 . but see IBN2002I . IIBi2003l . IN2003I and IM2003I for details. For us, the 
precise nature of this functional dependence does not matter. 

In fact, it is only in the computer science thesis IN2003I that the complete explicit 
parametric computation of the 8 families of differential forms uj, uj a , uj a , ip, tpg, (p a , ip a , 
ip together with the 8 invariant tensors S%%, i2| , T| 7 , Q a , L af) , P a(3 , H a and K a is 
achieved, in the case n = 2 and with the help of Maple. Although the task is really of 
impressive size, the author of this paper has the project of achieving manually the general 
computation for n > 2 variables (|M2005|). In IHa 19371 . the parametric computations 
are achieved completely only in the case n = 2. According to the author's experience, it 
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appears that the tensorial formalism helps to shorten importantly the size of the electronic 
computations and that the case n > 2 is not much more difficult than the case n = 2, as 
one may already observe by reading |M2004a|. This is why the project |M2005| seems 
to be an accessible task. A similar project would be to compute in a parametric way 
the structure equations obtained in |FeT9951 for the equivalence problem associated to a 
system of Newtonian particles; the computations are quite similar, indeed, as well as the 
computations of this paper as quite similar to the computations of | M2004a | . This will be 
achieved in the future, if time permits. 

At present, fortunately, there is a kind of "miracle": the interesting tensor Sg^ appears 
essentially at the beginning of the computation of the final {e}-structure. Thus, it is 
not necessary to go up to the end of the algorithm in order to deduce the final parametric 
expression of S%Z, and especially to characterize the maximally symmetric systems (1.2), 
i.e. those for which all the 8 tensors Sffi, R% , T^ 7 , L af3 , P a(3 , H a and K a vanish. 
In IM2003I . we obtained: 



It may be verified (cf. also IBi2003l ) that the general tensor is a certain "tensorial 
rotation" of its value "at the identity", namely its value when the matrix (1.10) is the 
identity; it suffices to put u := 1, :— 5p and u'^ := 5p in (1.15) above. Next, it may 
be verified that the vanishing of the value of at the identity (which is equivalent to 
the vanishing of the general So* ) provides a system of second linear second order partial 
differential equations involving only the partial derivatives Fy^y , ■ Finally, one verifies 
(" 16120 031) that this last system is equivalent to the fact that the F J1 J2 are of the specific 
cubic polynomial form (1.3). In conclusion, this provides a second (very summarized) 
proof of Theorem 1.1. 

In the remainder of the paper, we shall not speak anymore of the equivalence method 
and we will focus on describing precisely how to establish Theorem 1.1 using S. Lie's 
original techniques, as in IM2004al . 
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(LIFL, University of Lille 1), who implemented the Elie Cartan equivalence algorithm 
and who discovered the validity of Theorem 1.1 in the case n = 2, and also of Theo- 
rem 1.7 of IM2004al in the case m = 2. Without these discoveries, we would not have 
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1.18. Closing remark. After Theorem 1.1 was established, namely after the refer- 
ences IBN2002I . IN2003I . HBi2003l and IM2003I were completed, we discovered at the 
mathematical library of the Ecole Normale Superieure of Paris that Mohsen Hachtroudi, 
an Iranian student of Elie Cartan, also obtained a proof of Theorem 1.1 ( IHal937l . p. 53), 
sixty seven years ago. However, his computations follow the method of equivalence, 
in the spirit of his master, whereas we conduct ours in the spirit of Sophus Lie. Also, 
in IHal937l . M. Hachtroudi only refers to the abbreviated formulation (1.4) of the com- 
patibility conditions, and one does not find there the explicit and complete formulation of 
the systems (I'), (IT), (III') and (IV). Finally, the combinatorial formulas that we provide 
in (2.10), in (2.14) in (2.15), in (2.22), in (2.23), in (2.31) and in (2.35) below seem to b e 
new in the Lie theory of symmetries of partial differential equations (see also IM2004bl ). 



§2. Completely integrable systems of 
second order ordinary differential equations 

2.1. Prolongation of a point transformation to the second order jet space. Let K = R 

or C, let n£H, suppose n > 2, let x = (x , . . . , x n ) S K" and let y E K. According to 
the main assumption of Theorem 1.1, we have to consider a local K-analytic diffeomor- 
phism of the form 

(2.2) (a* , y) .— > (X^x^ , y), Y(x jl , y)) , 

which transforms the system (1.2) to the system Y X i lX i 2 = 0, 1 < ji, j% < n. Without 
loss of generality, we shall assume that this transformation is close to the identity. To 
obtain the precise expression (2.35) of the transformed system (1.2), we have to prolong 
the above diffeomorphism to the second order jet space. We introduce the coordinates 
{xi , y, y x ii , y x ii x i 2 ) on the second order jet space. Let 

d d <A F> 

(2-3) D k := -—r + y x k — + } y x k x i 



ox ay ^—^ ay. 



be the k-th total differentiation operator. According to |BK1989|, for the first order partial 
derivatives, one has the (implicit, compact) expression: 



(2.4) 





■I 








D x X n \ ( DiY 



D n X n / \ D n Y 



where (-)" 1 denotes the inverse matrix, which exists, since the transformation (2.2) 
is close to the identity. For the second order partial derivatives, again according 
to | BK1989 1, one has the (implicit, compact) expressions: 



(2.5) 





-( 







DxX 1 ■ 


■ D x X n \ 




( D 1 Y xi \ 


DnX 1 ■ 


■ D n X n J 




V D n Y XJ J 
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for 7 = 1, .... n. Let DX denote the matrix (DiX : '\\ < ?.'^ n , where i is the index of lines 

■' ' ' V 1 / l<i<n 

and j the index of columns, let Yx denote the column matrix (lx')i<i<n an< ^ l et DY be 
the column matrix (Z)jY) 1<i<n . 

By inspecting (2.5) above, we see that the equivalence between (i), (ii) and (iii) just 
below is obvious: 

Lemma 2.6. The following conditions are equivalent: 

(i) the differential equations Y X j x k — hold for 1 < j, k < n; 

(ii) the matrix equations D k (Yx) = hold for 1 < k < n; 

(iii) the matrix equations DX ■ D k (Yx) = hold for 1 < k < n; 

(iv) the matrix equations = D k {DX) ■ Yx — Dk(DY) hold for 1 < k < n. 

Formally, in the sequel, it will be more convenient to achieve the explicit computations 
starting from condition (iv), since no matrix inversion at all is involved in it. 



Proof. Indeed, applying the total differentiation operator Du to the matrix equation (2.4) 
written under the equivalent form = DX ■ Yx — DY , we get: 



(2.7) 



= D k (DX) -Yx+DX- D X {Y X ) - D k (DY), 



so that the equivalence between (iii) and (iv) is now clear. 



□ 



2.8. An explicit formula in the case n = 2. Thus, we can start to develope explicitely 
the matrix equations 



(2.9) 



= D k {DX) ■ Y x - D k (DY). 



In it, some huge formal expressions are hidden behind the symbol D k . Proceeding induc- 
tively, we start by examinating the case n = 2 thoroughly. By direct computations which 
require to be clever, we reconstitute some 3x3 determinants in the four (in fact three) 
developed equations (2.9). After some work, the first equation is: 



Vx 1 x 1 ■ 



(2.10) 



X x* X l* 

x\ x\ 



Y r 



X y 
Y, 



X l -X" 2 X 1 i 

x\ 1 

X L X L 



x\ 

X 1 



X x* X x^ 

x\ x\ 

X L X A 

Y x i Y X 2 

^x 1 ^x^x 1 

^-^1^1 



X lly 

x\ 

% y 



r, 



Y,. 



■< y 

% 

Y, 



X\ 

X A 
K2 



Y 1 

X 2 , x 

XX 



X x 2 

x\ 

Y^ 



X y 

Yy 
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Vx 1 Vx 1 



Vx 1 Vx 2 



Vx 1 Vx^ Vx 



Vx 1 Vxi Vx 




X x 2 

K 2 



x l 

X y 



This formula and the two next (2.22), (2.23) have been checked by Sylvain Neut and 
Michel Petitot with the help of Maple. 



2.11. Comparison with the coefficients of the second prolongation of a vector field. 

At present, it is useful to make an illuminating digression which will help us to devise 
what is the general form of the development of the equations (2.9). Consider an arbitrary 
vector field of the form 



(2.12) 



k=l 



dx k 



dy 



where the coefficients X k and y are functions of (x l ,y). According to IOL 19861 . 
IBK198 9 I, there exists a unique prolongation C^ 2 ' of this vector field to the second order 
jet space, of the form 



(2.13) 







ji=i 



Jl=l 32=1 



d 



where the coefficients Yj 1 , Yj 1 j 2 may be computed by means of some inductive devices 
explained in IOL 19861 . IBK1989I . In IGM2003I (see also ISu2001l for formulas which 
do not involve the Kronecker symbol), we obtained the following perfect formulas: 



71 

^01,32 = feiiJ! + Vx* 1 ' ~yn i2 V ^ ^ja ^x^y ~ ^xhxiz } 



fe 1= l 

n n 



(2.14) { 



X X y * ki Vx " 2 ' 

k 1 = l k 2 = l 



n n n 



Z Z Z Vx-i Vx"2 Vx-s 



fei=i fc 2 =i fe 3 =i 



for ji,j2 = 1 , • • ■ , n. The expression of Yj 1 does not matter for us here. Specifying this 
formula to the the case n = 2 and taking account of the symmetry Y1.2 = Y2.1 we get 
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the following three second order coefficients: 



(2.15) 



Yi,i = 



3V* 1 + Vx 1 ■ {2^i y - X*i xl } + y x 2 ■ {-X%i x i } + 
+ J/xi 2/xi • { y yy - 2 X^iy) + y x i y x 2 ■ {-2 X^ y } + 

+ Vx^ Vx^ Vx^ ■ {-Xyy} + Vx^ Vx 1 Vx 2 ' {~ X yy} , 

y x i x 2 + Vx 1 ■ {yx^y - x x i x 2} + y X 2 ■ {y x i y - x x i x2 } + 



Vx 2 



{y 



vv 



X X ly 



X x 2 y } + 



X 2 \ 
yyi 

X X 2 X 2 



+ Vx 1 Vx 1 ■ {-X^ y } +y x 

+ VX* VX* ■ {~X x ly} + 

+ Vx* Vx^ Vx 2 ■ {-Xy y } + y x i y x z Vx 2 ■ {- 

y x 2 x 2 + y x l ■ {~X^2 X 2} + Vx 2 ■ {2y x 2,y - X X 2 2 X 2} + 

+ Vxi Vx* ■ {-2 X^y} + y x 2 y x 2 ■ {y yy - 2 X 2 2y } + 

+ VX 1 Vx 2 VX 2 ■ {'Xyy} + Vx! Vx 2 VX* ■ {-Xyy} ■ 

We would like to mention that the computation of J2 ,1 < ji, _ 2, above is easier 
than the verification of (2.10). Based on the three formulas (2.15), we claim that we can 
guess the second and the third equations, which would be obtained by developing and 
by simplifying (2.9), namely with y x i x 2 and with y x 2 x 2 instead of y x i x 2 in (2.10). Our 
dictionary to translate from the first formula (2.15) to (2.10) may be described as follows. 
Begin with the Jacobian determinant 

X x l X x* X \ 
X\ X\ Xf, 



(2.16) 



Y, 













)■ 




1 






: 








of the change of coordinates (2.2). Since this change of coordinates is close to the identity, 
we may consider that the following Jacobian matrix approximation holds: 



(2.17) 



The jacobian matrix has three columns. There are six possible second order derivatives 
with respect to the variables (x 1 , x 2 , y), namely 

(2.18) {')x 1 x 1 ^ {')x 1 x 2 ^ {')x 2 x 2 : {')x 1 y: {'*)x 2 y: (')yy- 

In the Jacobian determinant (2.16), by replacing any one of the three columns of first 
order derivatives with a column of second order derivatives, we obtain exactly 3x6 = 18 
possible determinants. For instance, by replacing the third column by the second order 
derivative (-^ij, or the first column by the second order derivative (•) a a a .i, we get: 



(2.19) 



K 2 
K 2 

K 2 



Y 1 

Y 2 
A x^y 



Y, 



■ '■ L !l 



Y 1 

^x 1 x 1 
Y 2 

^x 1 x 1 



Y, 



X l 2 
X x 2 

Y X 2 



X l 

X 2 
Yu 



y 



We recover the two determinants appearing in the second line of (2.10). On the other 
hand, according to the approximation (2.17), these two determinants are essentially equal 
to 

i o xl ly X* 

1 x l v =Y x i v or to X 2 lxl 1 =X xlxl . 





(2.20) 



Y, 



x 1 y 



DIFFERENTIAL CHARACTERIZATION OF Y~ 



= 0, 1 < ji,j 2 < n > 1 



13 



Consequently, in the second line of (2.10), up to a change to calligraphic letters, we 
recover the coefficient 



(2.21) 



2 y x i y 

of y Xl in the expression of Yi,i in (2.15). In conclusion, we have discovered how to pass 
symbolically from the first equation (2.15) to the equation (2.10) and conversely. 

Translating the second equation (2.15), we deduce, without any further computation, 
that the second equation which would be obtained by developing (2.9) in length, is: 



(2.22) 




1 yy x y 

Using the third equation (2.15), we also deduce, without any further computation, that the 
third equation which would be obtained by developing (2.9) in length, is: 



(2.23) 



X x 2 

x l 2 

K.2 



Y X 2 X 2 

x l 



X v 
X y 
Yy 

K 

x l 

Y x < 

X x 2 
X x 2 

Y»2 



x l 




yl 

^ X 2 X 2 


x l 


K 2 


Y 2 

^ X 2 X 2 


Y x i 


Y X 2 


Y X 2 X 2 



X y 
Yy 

^ X 2 y 

Y 2 
^ X 2 y 



Y, 



x - y 



X x^ 

x l 



X 

Yy 
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x\ 

x y 



X l 



Y 



Y, 



V 



Vx 1 Vx 2 



Vx 2 Vx 2 



Vx 1 Vx 2 Vx 



Vx 2 Vx 2 Vx 2 



2.24. Appropriate formalism. To describe the combinatorics underlying formu- 
las (2.10), (2.22) and (2.23), as in |M2004a|, let us introduce the following notation for 
the Jacobian determinant: 

xl xl> xl 

v2 v2 v2 
A_i A„ 2 A... 



(2.25) 



A(*Vlz/) 



Y,. 



Here, in the notation A(x 1 |a; 2 |y), the three spaces between the two vertical lines | refer 
to the three columns of the Jacobian determinant, and the terms x 1 , x 2 , y in (a; 1 |x 2 |j/) 
designate the partial derivatives appearing in each column. Accordingly, in the following 
two examples of modified Jacobian determinants: 



(2.26) 



A(x 1 x 2 \x 2 \v) 



A(x 1 |x 2 k 1 y) 





Y 1 

^x 1 x 2 


x l 


x\ 




Y 2 

A x^x 2 


X X 2 


X y 




Yx 1 x 2 


Yx 2 


Yy 




xl 


X X 2 


X ly 






x % 






Y x i 


Y X 2 


Y i 

A x L y 



and 



we simply mean which column of first order derivatives is replaced by a column of second 
order derivatives in the original Jacobian determinant. 

As there are 6 possible second order derivatives (-^i^i, (-) x i x 2, (-) x i x y, (-) X 2 X 2, (-) x 2 y 
and (-)yy together with 3 columns, we obtain 3 x 6 = 18 possible modified Jacobian 
determinants: 



(2.27) 



A(x 1 x 1 \x 2 \y) 

A(x x x 2 \x 2 \y) 

A{x l y\x 2 \y) 

A(x 2 x 2 \x 2 \y) 

A(x 2 y\x 2 \y) 

Myy\x 2 \y) 



A(x L \x 1 x 1 \y) 
A(o; 1 |a; 1 a; 2 |y) 

A(x 1 \x 2 x 2 \y) 
A(x 1 \x 2 y\y) 



A(a; 1 |a; 2 
Aix^x 2 
Aix^x 2 
A(x 1 \x 2 
A(xV 



\x x x x ) 
\x l x 2 ) 

\x x y) 

\X X ) 

\x 2 y) 
\yy)- 



Next, we observe that if we want to solve with respect to y x i x i in (2.10), with respect 
to y x i x 2 in (2.22) and with respect to y x 2 x 2 in (2.23), we have to divide by the Jacobian 
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determinant A(x 1 \x 2 \y). Consequently, we introduce 18 new square functions as follows: 
(2.28) 





A^xVly) 
A (a; 1 a; 2 |j/) 


Q^x 2 


A(a; 1 |a: 2 |y) 




A(a; 1 j/|a; 2 |2/) 
' A(x!|a; 2 |y) 




A(x 2 arV|y) 
A(s 1 a; 2 |y) 




A(aA/|:E 2 |i/) 
' A(a: 1 |a; 2 |y) 


D 1 
'— 'vv 


A(yj/|a; 2 |2/) 
■ A(x 1 |a; 2 |y) 




A{x 1 \x 2 \y) 


d^x 2 


A(a: 1 |a; 1 a; 2 |y) 
A{x 1 \x 2 \y) 




A(a: 1 |a: 1 y|y) 
' A(x 1 |i 2 |y) 




A(a; 1 \x 2 \y) 


°* 2 y 


A(x 1 \x 2 y\y) 
A(a- 1 |a; 2 |j/) 


□ 2 


• A(x^\x 2 \y) 




A(a; 1 |x 2 |a; 1 a; 1 ) 
A(x 1 \x 2 \y) 


D^X 2 


A(a; 1 |a; 2 |a: 1 a; 2 ) 
A(x 1 |x 2 |y) 




A(a; 1 |a; 2 |a; 1 y) 
A(a; 1 |a; 2 |y) 


D^x 2 

V 


A(x 1 \x 2 \y) 


□x 2 » 


A(x 1 \x 2 \x 2 y) 
A(x 1 \x 2 \y) 


□ 3 
'—'vv 


Aix 1 \x 2 \yy) 
' A(x!|a; 2 |y) ' 



Thanks to these notations, we can rewrite the three equations (2.10), (2.22) and (2.23) 
in a more compact style. 

Lemma 2.29. A completely integmble system of three second order partial differential 
equations 

!Vx^xA x ) = F 1 ' 1 {x 1 ,x 2 ,y{x),y x i{x),y x 2{x)) , 
Vxwix) = F 1 ' 2 ip 1 ,x 2 ,y{x),y x i{x) 1 y x 2{x)) , 
Vx 2 x 2 (x) = F 2 ' 2 (x 1 ,x 2 ,y(x),y x i(x),y x 2(x)) , 

is equivalent to the simplest system Yx x x x = 0> Y x i X 2 = 0, Yx 2 X 2 = 0, if an d only 
if there exist local ¥L-analy tic functions X , X 2 , Y such that it may be written under the 
specific form: 

' Vx^x 1 = -QxV + Vx 1 ■ (-2 0^1^ + O x i x i) + y x 2 ■ (dli x i) + 
+ Vx* Vxi ■ (~ a 'yy + 2 G *i y ) + y x i y x 2 ■ (2 D 2 i y ) + 
+ Vx 1 Vx 1 Vx 1 ■ ( a yy) + Vx^ Vx 1 Vx 2 ■ iply) , 

Hx^x^ = -n^x 2 + Vx± ■ {-^liy + n^ix 2 ) + Vx 2 ■ + o 2 i a ,2) + 

+ 2/xl • ( D i 2 y) + 2^ ^ 2 ' {- U lv + G ily + G * 2 y) + 

+ Vx* Vx^ ■ (n^iy) + Vx 1 Vx 1 Vx 2 ■ + 2/x 1 2/x 2 Vx? ■ , 

J/A 2 = -d^ 2 + fci ■ ( n i 2 z 2 ) + 2^2 ■ (-2 U x 2 y + 2 X 2 X 2) + 
+ Vx* Vx 2 ■ (2 Oily) + Vx 2 Vx* ■ (~ a yy + 2 ^ 2 y ) + 
+ Vx 1 Vx? Vx 2 ■ (n^y) + Vx 2 Vx* Vx 2 ■ (Oyy) ■ 

2.32. General formulas. The formal dictionary between the original determinantial for- 
mulas (2.10), (2.22), (2.23), between the coefficients (2.15) of the second order prolon- 
gation of a vector field and between the new square formulas (2.31) above is evident. 
Consequently, without any computation, just by translating the family of formulas (2.14), 
we may deduce the exact formulation of the desired generalization of Lemma 2.29 above. 

Lemma 2.33. A completely integrable system of second order partial differential equa- 
tions of the form 

(2.34) y x iixJ2{x) = F n ' j2 (x,y(x),y x i(x), . . .,y x «(x)) , j±, j 2 = l,...n, 



(2.31) 
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is equivalent to the simplest system Y X ii x?2 = 0, ji, j% = 1, • • . , n, if and only if there 
exist local ^.-analytic functions X , Y such that it may be written under the specific form: 



(2.35) 

Vx^x^ 



+ V y xkl • ( ( U\ 39 - S kl □"t 1 - 6 kl □"jt 1 



fei=i 
k- 1 



LT 1 

x'2y 



- <5 Kl LT 
2 



Vxii 



□ fc i _ I ,J fe i □«+! 

ins 2 J 1 »f 



+y x31 y xi2 ■ (□*) } 



Of course, to define the square functions in the context of n > 2 independent variables 
(or, x 2 , . . . , x n ), we introduce the Jacobian determinant 



(2.36) 



A(x 1 \x 2 \---\x n \y) :-- 



XL. 



v n Vfi 



\y), 



together with its modifications 

(2.37) A(x 1 \---\ kl x jl x j2 \ 

in which the fci-th column of partial first order derivatives I* 1 x kl | is replaced by the col- 
umn | fel x^x^l of partial derivatives. Here, the indices fci, juj2 satisfy 1 < fci,ji,j2 < 
n + 1, with the convention that we adopt the notational equivalence 



(2.38) 



n+l 



y 



This convention will be convenient to write some of our general formulas in the sequel. 

As we promised to only summarize the proof of Theorem 1 . 1 in this paper, we will not 
reproduce the proof of Lemma 2.37 from IIM2U03I . Involving linear algebra considera- 
tions, this proof is similar to (and in fact slightly simpler than) the proof of the analogous 
Lemma 3.32 in IM2004al . 



§3. First and second auxiliary system 

3.1. Functions G 7 , H4 1 „• , L kl and M kl . To discover the four families of functions 
appearing in the statement of Theorem 1.1, by comparing (2.35) and (1.3), it suffices (of 
course) to set: 



(3.2) 



~ji ,32 

tiki 

31 ,32 

L kl 

31 

M kl 



■= -LT+ 1 ■ 

xHxJ2> 

:= n\ - S kl LT^t 1 - 5 kl LT+ 1 , 

X 31 X 32 Jl Xl2y J2 x^ly' 
o ri^i rfci r— in+1 

:= D k i. 



Consequently, we have shown the "only if" part of Theorem 1.1, which is the easiest 
implication. 

To establish the "if" part, by far the most difficult implication, the very main lemma 
can be stated as follows. 
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Lemma 3.3. The partial differerential relations (I'), (II'), (HI') an^ (IV') which express 
in length the compatibility conditions for the system (1.3) are necessary and sufficient for 
the existence of functions X, Y of {x l1 , y) satisfying the second order nonlinear system 
of partial differential equations (3.2) above. 

Indeed, the collection of equations (3.2) is a system of partial differential equations 
with unknowns X 1 , Y, by virtue of the definition of the square functions. 



3.4. First auxiliary system. To proceed further, we observe that there are (m + 1) more 
square functions than functions Gj 1} j 2 , H-^^, and M kl . Indeed, a simple counting 
yields: 



(3.5) 



whereas 



(3.6) 



#{G J1J2 } = 



n 2 (n+ 1) 



n(n + 1) 



+1 n(n + l) 

= i, 

#i H h,hI - — 2 — ' 



#{M fcl } = n. 



Here, the indices j\, ji, k\ satisfy 1 < jx,j2,ki < n. Similarly as in |M2004a|, to 
transform the system (3.2) in a true complete system, let us introduce functions Jo of 
(a;' 1 , y), where 1 < j\ ,j2,k\ < n + 1, which satisfy the symmetry IT^ 1 - a 
let us introduce the following first auxiliary system: 




LT+ 1 

x J1 y 



fei 

ji,n+l> 

n+l 



n fei = n fci 

1 — l yy AJ -n+l,n+l' 

n«+ 1 — TT"+ 1 

LJ TO — 11 ri+l,7i+l- 



(3.7) 

1 ji,n+l> 

It is complete. The necessary and sufficient conditions for the existence of solutions X 1 , 
Y follow by cross differentiations. 



Lemma 3.8. For all ji, j'2, J3, fci 

relations 

(3.9) 

{ D x^x") x j 3 ~ (P^ 1 ^ 3 ) x i2 = 



1,2,.. 



1, we have the cross differentiation 



n+l 

£ 

fc 2 =i 



n+l 



□ 



□ 



fci 



Z ^ X J lx 3 3 X 3 2x* 



The proof of this lemma is exactly the same as the proof of Lemma 3.40 in IM2004al . 

As a direct consequence, we deduce that a necessary and sufficient condition for the 
existence of solutions IT^ 1 - a to the first auxiliary system is that they satisfy the following 
compatibility partial differential relations: 



(3.10) 



an 



fci 



an 



JU3 



U kl , 

33M 



for all jx,h,h, fa = 1, ... , n + 1. 

We shall have to specify this system in length according to the splitting {1,2, ... ,n} 
and {n + 1} of the indices of coordinates. We obtain six families of equations equivalent 
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to (3.10) just above: 
(3.11) 



\ h,h) x 33 \ jl ,33 1 x 32 Z_^l 31,32 hM 31.32 33,1+1" 



k 2 = l 



1 \ v y\k 2 TT™+ 1 J- TT n + 1 TT n + 1 

fc 2 = l 

n 

(H^jUy ~~ (-^illn+l)^ = ~ X] n il 2 j2 n n+l,fc 2 ~ n ilj2 n "+l.n + l + 

fc 2 =l 

n 

1 \ v y\k 2 TT™+ 1 _L TT n + 1 TT n + 1 
i " u ji,n+l u jj,iij " l " il 2i,n+l iX 33,n+l) 



k 2 = l 



(T\ n+1 ) - (TT™ +1 ) — — TT fc2 TT™ +1 - TT" +1 TT n+1 4- 

( lL j 1 ,n+l) y \ ll n+l,n+l) x 3 1 — ^ 11 j 1 ,n+l LL n+l,k 2 AA ji,n+l 1L n + l.n+l ^ 

k 2 = l ^ 
n 

1 \ y\k 2 TT n + 1 -I- TT™+ 1 TT"+ 1 

fe 2 =l 



TT fel 1 -I TT fel I — - TT fc2 TT fel - TT" +1 TT fcl 4- 
"./. ,32 ) ja \ U j ,h ) xi2 ~ 2^ ,h LL 33 M L1 ji ,32 U js ,n+l + 

ka=l 

n 

1 yik 2 i-rfci , TTn+1 rrfei 

+ 2^ li ii,i 3 li j 2 ,fe 2 +AA ii,j3 AA j2,n+i' 
fc 2 =i 

n 

TT fcl \ — I T\ kl \ — - TT fc2 TT fel — TT n+1 TT fel 4 
M'Jv \ jl ' n+1 )xl2 ~ ii ii,J 2 AA n+l,fe 2 A1 n,j2 AA n+l,n+l + 

fc 2 =l 



fci 

j' 2 ,n+l' 

fc 2 =l 



4- TT fe2 TT fel 4- TT n+1 TT 

fc 2 =l 

n 

( n t + l) -( n " +W ) Kjl =-E U tn + lK 

y 7, 1 



fci _ rr n + 1 TT^ 1 
+l,fc 2 AA ji,n+l LL n+l,n+l 

k 2 = l 



En' 2 TT fcl 4 TT n+1 TT' £i 

11 n+l, ) i+l 1J -ii,fc2 ' AA n+l,n+l LL ji,n+V 

k 2 = l 



where the indices ji, j 2 , J3, fci vary in the set {1, 2, 1, . . . , n}. 



3.12. Principal unknowns. As there are (m + 1) more square (or Pi) functions than the 
functions Gj lt j 2 , - 9 , L - 1 and M fel , we cannot invert directly the linear system (3.2). 
To quasi-inverse it, we choose the (m + 1) specific square functions 



(3.13) 



Q :—O x i xl , 9 2 O x2x2 , 
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calling them principal unknowns, and we get the quasi-inversion: 
(3.14) 

Jjfei jyAj _ 1 gkj, H J2 jfcl ^1 I I ^1 0J2 , I jfel @J1 

Ji J2 zJiz^ 11 n ,] 2 2 Jl ■ j2 ' 3a 2 32 3l ' 3l 2 Jl 2 32 ' 

-rrn+1 _ rnn+1 _ _r~i 
TT" +1 — n " + 1 — — - H n 4- i B- 71 



3.15. Second auxiliary system. Replacing the five families of functions IF^ 1 - 2 , IT?* n+v 



H-n+i ,n+i> n ii + / 2 ' n ii + n+i ^ their values obtained in (3.14) just above together with the 
principal unknowns 



(3.16) 



rpi — p>h TT™ +1 — B™ +1 

Ji.ii ~~ ' 11 n+l,n+l — " > 



in the six equations (3.11)i, (3.11) 2 , (3.11) 3 , (3.11) 4 , (3.11) 5 and (3.11) 6 , after hard 
computations that we will not reproduce here, we obtain six families of equations. From 
now on, we abbreviate every sum as Y^k • 

Firstly: 

(3.17) = G jld2tX j 3 - G juh ^ j2 + ^2 G h M H ji,j 2 - E G i 2 > fc i H h,h • 



This is (F) of Theorem 1.1. Just above and below, we plainly underline the monomials 
involving a first order derivative. Secondly: 



( e^i 

xJ2 



E G nM L t + \ H ln H ln E H n-» H tW 



(3.18) ^ 



fci 



fci 



Q. . B ,l+1 H jl B 32 H J2 B J1 + H kl B fcl - 



ki 



Thirdly: 



(3.19) 



_ J1 Q J2 



_e n+1 + - 6 J1 = - H n - - 

a" ^ 2 w 2 Ji' 3 !'" 



E G ;^ M * 1 + iE<* 1 4 1 + 

- fff . 6 ,l+1 - - V L^ 1 B fel - - B jl B' 1+1 . 

4 JUi 4 Ji 4 
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Fourtly: 
(3.20) 

' I s ki e h - - s ki e h + - 5 ki e ji - - s ki e Ji = 

= -# fcl ■ +H kl . - - c5 fel iJ i3 ■ + - (5 fcl # J2 

^ //' + ls kl H jl . ,- 
2 ji,ji,x« 2 j2 JiiJii"^ 3 



^ r fei 4-^r* r fei ^ A*-' 1 ft? 1 Pf j2 _i_ \ A fci ft-? 1 ft-? 3 

Errk-2 rrfci , \ rrfta rrfci _ 1 rfei rrb rrfca , 1 rrfc 2 rrfe 

fl 3U2 n jj,kj + ^ii J3 n jj,kj 2 J2 ••' fe 2,fc2 + 2 « j2 fc 2,fc 2 

&2 &2 

- \ s n H tn ® js • K *Li (_) K ^,,3 © J • t< ( - ),: - 

_ I s ki V ij fe2 e fc2 + - 6 ki V h* 2 e* 2 - 

fcl fel 

- e n & n + -s ki e 31 e h . 

Fifthly: 
(3.21) 

' I A fe i + 1 A fe i e^i _ I s ki e n+1 = 

2 ji v T 2 J2 y 2 ■ Jl K 



G 71 *, M fel + -V L* 2 - - V L? 1 - - 6 kl Y H k \ L 

k 2 ^2 k-2 

- - s ki h] 1 . e ,l+1 + - s ki V L k2 e fe2 + - s ki e 31 e n+1 . 

A 32 3l,3l A 32 / i 3l A 32 



Sixthly: 
(3.22) 



+ 2 V H kl , M fc2 - V H k2 , M l! --V L^ 2 Li 1 - 

Al2 Al2 fi2 

+ s ki y M k2 e fe2 + - s ki e n+1 e ,l+1 

jl / .j y 31 

k 2 



3.23. Solving 6^ 2 , G^ 1 , O™^ 1 and 6™ +1 . From the six families of equations (3.17), 

(3.18), (3.19), (3.20), (3.21) and (3.22), we can solve Q j x ) 2 , Q^ 1 , O™^ 1 and 9™ +1 . Not 

mentioning the (hard) intermediate computations, we obtain firstly: 

(3.24) 

^2 = ~ 2 G 3i,32,y + H jl n , x i2 + X! G h,l L ji + 2 H jl,3i H il,32 - Yl H n,32 H l,l' 

I I 



Secondly: 
(3.25) 

+ 1 E ^ - 1 E H Ln V - 1 K,, + ^E4, 

( Z z 

+ I e ji e«+i. 

Thirdly: 
(3.26) 



3 ^ jwi * 2 2 
z z 



1 ji Q n+1 



Fourtly: 
(3.27) 



2 

z z z 



y m ! e' + - e n+1 e n+1 



2 

These four families of partial differential equations constitute the second auxiliary system. 
By replacing these solutions in the three remaining families of equations (3.20), (3.21) 
and (3.22), we obtain supplementary equations (which we do not copy) that are direct 
consequences of (F), (IF), (III'), (IV). 

To complete the proof of the main Lemma 3.3 above, it suffices now to establish the 
first implication of the following list, since the other three have been already established. 

• Some given functions Gj 1 j 2 , H^j, and M kl of (x ll ,y) satisfy the four 
families of partial differential equations (F ), (IF ), (IIF ) and (IV ) of Theorem 1.1. 

• There exist functions jl , Q n+1 satisfying the second auxiliary system (3.24), 
(3.25), (3.26) and (3.27). 

• These solution functions J1 , 8™ +1 satisfy the six families of partial differential 
equations (3.17), (3.18), (3.19), (3.20), (3.21) and (3.22). 

• There exist functions 11** - 2 of (x ll ,y), 1 < ji,j2,ki < m + 1, satisfying the 
first auxiliary system (3.7) of partial differential equations. 

• There exist functions X 1 ' 2 , Y of (x ll ,y) transforming the system y x ji x i2 = 
F jl ' j2 (x h ,y, y x i 2 ), ji,j2 = 1, to the simplest system Y X j 1X j 2 = 0, 
ji,ji = l,...,n. 
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3.28. Compatibility conditions for the second auxiliary system. We notice that the 
second auxiliary system is also a complete system. Thus, to establish the first above 
implication, it suffices to show that the four families of compatibility conditions: 



k o=(oS 1 )»-(eS +1 U. 



are a consequence of (F), (I"), (III'), (IV). 

For instance, in (3.29)i, replacing 9^ 2 by its expression (3.24), differentiating it with 
respect to x 33 , replacing J ^ 3 by its expression (3.24), differentiating it with respect to 
x 32 and substracting, we get: 
(3.30) 

+ l <") ' , + l <-)' H'\ - I (-)' , O ' 1 - i (-)' <-) ' •.. 



- - // ' . (-)•' - - //; , (-)•' + - //' . h (-)• + - //; , Q j % 

- - // ' . „ e* - - //; , (-)•' + - //•'• . ia e* + - //•' , e 31 J2 

2 J2,J2,K 3 3 2 J2.J2 X J 3 2 33,33, x 2 2 J3j3 X J 2 



y a;32 

z i Z Z 

+ I // ■•: . // •' ■ . + I //•' . // " . _ I // ■' . // •" . _ I // ■' . // " 

2 ii.ii^ 33 J2j2 2 31,31 J2,i2,a; J 3 2 jiji,^ 32 J3,j3 2 33,33^2 



Y H n-J2-x^ "ij - Y 11 [ ... #U,^3 + J! H jnJ3^ 2 ' + Y "'■ H U 



X 3 2 



Y G h,l,x i3 L h + Y l ' " J L jl,x j 3 ~ Y G 33,l,xi2 l'r - Y (,, - J L . 



33,1 lj j 1 ,x^2 ■ 



Next, replacing the twelve first order partial derivatives underlined just above: 



(3.31) 



Q J1 

ZJ3 ■ 

Z->3 ■ 



& 32 

^J3 ■ 

^32 ■■ 



Q 31 

gJ2 ' 

Qn+1 

X^3 ' 



& 33 

gJ2 ' 

Qn+1 

' 



& 32 

X^3 ' 



& 33 

x^2 : 



e 



X33 ' 



e 



Xl2 ■ 
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by their values issued from (3.24), (3.26) and adapting the summation indices, we get the 

explicit developed form of the first family of compatibility conditions (3.29) 1: 

(3.32) 

=? = -2Gj lt j 2 ,x*3y +^ G 3u j 3 ,xi2y- 



~ Y, •< ■•• L h + Y, Gj 2t i, x i 3 L l n Gj lt j 2<y II fl,. + G n .j 3 .y Hf^j, 

_t I 

- 2 Y ■ H k,h ■ 2 Y "' ■ - Y " J ■ H U + Y H Lh,*™ "' ' 

; i i i 

2 2 2 2 

_J [ 

2 ; 2 4 >— ^ i 

~ 3 Cji j2 Cj3 J3 M 3 + - Gj u j 3 (•',,.,, M n - G n j 2 Gj 3> t M + 

4 



Y G h>h ( 'J-' }l ' 2 Y ( "j. j L h + 2E (: >' J H n-h L h' 
l 1 1 

o Y (! ' J H h,h L h + \Y (: > j H h,h L h ~> Y G ii>i* L k+ 
1 1 1 

' 2 Y G iuh H l,l L h 3 Y G n-n Hi' • L j 3 + 3E G h,ia H j 2 ,i L h~ 
1 1 1 

~ 3 G h,js 11 L T + 3 G h,h 11 [..j, L l 3 ~ 

Y Y f '"'•<" H h,h L i ' Y Y ( ' ' H k,h L 'i 
I p I p 

-YY H Lh H lh h p,p + YY H Lh H lj2 h p,p- 
1 p 1 p 

Lemma 3.33. (|M2003|, |M2004a|) This first family of compatibility conditions for the 
second auxiliary system obtained by developing (3.29)i in length, together with the three 
remaining families obtained by developing (3.29)2, (3.29)3, (3.29)4 in length, are con- 
sequences, by linear combinations and by differentiations, of (V), (IV), (III'), (IV), of 
Theorem 1.1. 

The proof of Theorem 1 . 1 is complete. □ 
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